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Abstract 

Botanical primary metabolites extensively exist in herbal medicine injections (HMIs), but often were ignored to control. With 
the limitation of bias towards hydrophilic substances, the primary metabolites with strong polarity, such as saccharides, 
amino acids and organic acids, are usually difficult to detect by the routinely applied reversed-phase chromatographic 
fingerprint technology. In this study, a proton nuclear magnetic resonance {^H NMR) profiling method was developed for 
efficient identification and quantification of small polar molecules, mostly primary metabolites in HMIs. A commonly used 
medicine, Danhong injection (DHI), was employed as a model. With the developed method, 23 primary metabolites 
together with 7 polyphenolic acids were simultaneously identified, of which 13 metabolites with fully separated proton 
signals were quantified and employed for further multivariate quality control assay. The quantitative NMR method was 
validated with good linearity, precision, repeatability, stability and accuracy. Based on independence principal component 
analysis (IPCA), the contents of 13 metabolites were characterized and dimensionally reduced into the first two 
independence principal components (IPCs). IPCl and IPC2 were then used to calculate the upper control limits (with 99% 
confidence ellipsoids) of and Hotelling control charts. Through the constructed upper control limits, the proposed 
method was successfully applied to 36 batches of DHI to examine the out-of control sample with the perturbed levels of 
succinate, malonate, glucose, fructose, salvianic acid and protocatechuic aldehyde. The integrated strategy has provided a 
reliable approach to identify and quantify multiple polar metabolites of DHI in one fingerprinting spectrum, and it has also 
assisted in the establishment of IPCA models for the multivariate statistical evaluation of HMIs. 
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Introduction 

Metabolic profiling is essential to ensure the quality, consisten- 
cy, safety and efficacy of herbal medicine products; especially for 
the injection dosage form. Water extraction or decoction is the 
most favored method of preparation of herbal medicine injections 
(HMIs). Saccharides, amino acids, organic acids, and other 
primary metabolites are unavoidably extracted along with targeted 
secondary metabolites during the process of HMIs, such as 
QingkaUing injection [1], Danshen injection [2], Guanxinning 
injection [3], and Shuxuetong injection [4]. In general, amino 
acids are considered to provide tonic activities and act as key 
regulators of nutrient metabolism, and polysaccharides are 
believed to be one of the most important constitutes in some 
herbal materials for pharmacological activities. Some monosac- 
charides have the suppressive effect on cell-mediated immune 



reactions [5]. However, these primary metabolites in HMIs are 
often ignored to detect and set corresponding quality criteria in 
China Pharmacopeia and national standard. Due to the strong 
polarity and hydrophUicity, conventional reversed-phase perfor- 
mance liquid chromatography (HPLC) fingerprint [6,7]and other 
analytical methodologies, including anion exchange chromatog- 
raphy [8], gas chromatography [9,10], and capillary electropho- 
resis [11,12], are limited to separate and detect primary 
metabolites unless with complicated pretreatment, derivatization 
reagents and laborious preparation procedures. To delineate 
various class metabolites in HMIs, different types of fingerprints 
are necessary for a holistic quality evaluation, which are difficult to 
realize during the practical industry processing. Therefore, a 
simple and fast approach is required to be capable of detecting 
saccharides, amino acids, and organic acids, together with mainly 
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bioactive secondary metabolites in one fingerprinting spectrum 
simultaneously. 

Proton nuclear magnetic resonance ( HNMR) spectroscopy 
provides access to detect all proton-bearing compounds, almost 
irrespective of the chemical compound class. Because of the signal 
intensity directly proportional to the number of nucleus contrib- 
uting to a specific resonance [13], 'H NMR method achieved the 
identification and quantification of metabolites in a one-step 
acquisition. With simple sample preparation, H NMR facilitates 
high-throughput analysis for metabolic studies and quality control 
of various food [14-17] and herbal materials [18-20]. As a 
universal technique with the simplicity and rapidity of implemen- 
tation, ' H NMR has expansive prospect of application in profiling 
of polar metabolites of HMIs. 

Since the complexity of chemical composition, HMIs cannot be 
completely represented by a limited number of certain bioactive 
compounds [21]. To extract the feature variables, principal 
component analysis (PCA) [22] and independent component 
analysis (ICA) [23] are classical tools to reduce the dimension of 
multivariate. The components (PCs) in PCA method are mutually 
orthogonal, while ICA method contains the components to be 
statistically independent [24]. ICA has been found to be a 
successful alternative to PCA in eliminating the overlapping 
information between the components [25]. However, ICA faces 
some limitations due to some instability [26], the choice of number 
of components to extract and high dimensionality [27]. As a 
consequence, independent principal component analysis (IPCA) 
was proposed by Yao et al. in 2012 [28] to use PCA as a pre- 
processing step to reduce the dimension of the data, and then use 
ICA as a denoising process of PCA to separate relevant 
information. On simulation studies and real data sets, IPCA 



offered a better visualization of the data than ICA and with a 
smaller number of components than PCA. Owing to the benefit to 
generate denoised the loading vectors, we attempted to employ 
IPCA method to construct and HoteUiiig T^ control charts [29] 
for multivariate statistical analysis. 

Danhong injection (DHI) is a patent injection made from the 
extracts of Radix Salviae Miltiorrhizae and Flos Cmihami [30] . It 
has been widely used for the prevention and treatment of 
cardiovascular and cerebrovascular diseases in clinic [31-33]. In 
our previous work, ultra-performance liquid chromatography 
(UPLC) coupled with UV detection was adopted to identify 
1 1 polyphenolic acids in DHI [34]. However, the total weight of 
identified constituents accounted for only a low proportion (about 
10%) of the solid content in DHI. 

In this study, we describe a strategy to detect more hydrophilic 
primary metabolites in DHI based on quantitative 'H NMR 
spectroscopy. The absolute concentration of identified metabolites 
was calculated by using internal standard method. Linearity, 
precision, repeatability, stability and accuracy were carried out to 
validate the method. The contents of polar metabolites were 
further evaluated by the multivariate analysis tool IPCA to 
establish and HoteUing T^ control charts. 

Materials and Methods 

Materials and Chemicals 

Thirty-six batches of DHI manufactured in 2011, 2012 and 
2013 were provided by Heze Buchang Pharmaceutical Co. Ltd 
(Heze, China). The standards of valine, threonine, alanine, 
pyroglutamate, procatechuic aldehyde and asparagine were 
purchased from the National Institute for Food and Drug Control 
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Figure 1. Representative NiVlR spectra of Danhong injection. Peaks: 1, Isoleucine; 2, Leucine 3, Valine; 4, Threonine; 5, Alanine; 6, Acetate; 
7, Proline; 8, Pyroglutamate; 9, Succinate; 10, Asparagine; 11, Malonate; 12, Glucose; 13, Galactose; 14, Arabinose; 15, Fructose 16, Rhamnose 17, 
Rutinose 1 8, Rutinulose; 1 9, Salvianic acid; 20, Salvianolic acid B; 21 , Rosmarinic acid; 22, Lithospermic acid 23, Salvianolic acid A; 24, Procatechuic acid; 
25, Procatechuic aldehyde; 26, 4-Hydroxybenzoic acid; 27, 4-Hydroxycinnamic acid; 28, Uridine; 29, Formate; 30, 5-(Hydroxymethyl)-2-furaldehyde. 
doi:1 0.1 371/journal.pone.01 0541 2.g001 
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Table 1. Chemical shifts for 30 metabolites identified in the NMR spectra of Danhong injection. 



Peak 


Metabolites 


Assignment 


ifH (multiplicity") 


Assigned with 


1 


Isoleucine 


S-CH3, y'-CHg, y-CH2, ^-CH, a-CH 


0.96 (t, 6.1 Hz), 1.01 (d, 6.8 Hz), 1.26 (rr), 1.47 (m), 
1 .96 (mj, D.Ob (m) 


COSY, TOCSY, J-res, 
HSQC 


2 


Leucine 


^ S' f~LJ 0 /"LI .. /"LI 

0, 0 -LH3, />-Ln2,; -Ln 


u.yD (mj, 1.0/ (m), D./2 (m) 


ludY, IuldY, j-res, 
HSQC 


3 


Valine 


y -CH3, j'-CHs, p-LH, a-LH 


1.00 {d, 7.0 Hz), 1.05 (d, 7.0 Hz), 2.26 (m), 3.60 (d, 4.2 Hz) 


/"/^C\/ T/^/"C\/ / 

COSY, TOCSY, J-res, 
HSQC 


4 


Threonine 


/"II /? /"1 1 /"II 
T-CH3, p-CH, K-CH 


1.36 {d, 6.5 Hz), 4.25 (m), 3.58 (d, 6.4 Hz) 


COSY, TOLbY, J-res, 
HSQC 


5 


Alanine 


^-CHj, 2-CH 


1.48 {d, 7.5 Hz), 3.78 (m) 


COSY, TOCSY, J-res, 
HSQC 


6 


Acetate 


CH3 


1.93 {s} 


J-res, HSQC 


7 


Proline 


y-CHj, 15-CH2, J9-CH2, a-CH 


2.00 (m), 3.34 (m), 3.42 (m), 2.35 (m),4.14 (m) 


COSY, TOCSY, J-res, 
HSQC 


8 


Pyroglutamate 


y-CHj, /i-CHj, s(-CH 


2.04 (m), 2.51 (m), 2.41 (m), 4.18 (m) 


COSY, TOCSY, J-res, 
HSQC 


9 


Succinate 


CH2 


2.44 (s) 


J-res, HSQC 


1 0 


Asparagine 


/?-CH2, a-CH 


2./6 (uu, 14.0, 0.2 Hz), 2.9/ lad, 14.U, 4.3 Hzj, 4.0U (m) 


COdY, IOCdY, J-res, 
HSQC 


1 1 


Malonate 


CHj 


3.19 (s) 


J-res, HbQC 


12 


Glucose 


aCiH, aCjH, aCjH, 11C4H, aCjH, 

O/" LI (3/~ LI A?/" LI 

pt-iH, /5C2n, pCjH 


5.24 {d, 3.6 Hz), 3.55 (m), 3.72 (m), 3.42 (m), 3.85 (m), 
4.DD tu, /.o Hz), 3.25 (dd, /.o, 9.2 Hz), 3.48 (mj 


COSY, TOCSY, J-res, 
HSQC 


13 


Galactose 


aCiH, aCiH, aCsH, ^CiH, /iCjH, 

O/" 1 1 O/" 1 1 

PC3H, PC4H 


5.27 (d, 3.8 Hz), 3.83 (m), 4.02 (m), 4.59 (d, 7.8 Hz), 3.50 (m), 
3.65 (m), 3.93 (m) 


COSY, TOCSY, J-res, 
HSQC 


14 


Arabinose 


aCjH, cfCjH, aCjH, /JC,H, /iCjH, 
/JCsH 


5.25 {d, 3.6 Hz), 3.82 (m), 4.02 (m), 4.52 {d, 7.8 Hz), 3.51 (m), 
3.68 (m) 


COSY, TOCSY, J-res, 
HSQC 


15 


Fructose 


CiH, C3H, C4H, C5H, QH 


3.58 (m), 3.69 (m), 3.82 (m), 4.11 (m), 3.90 (m), 4.11 (m), 3.82 (m), 
4.00 (m), 3.69 (m), 3.82 (m), 4.02 (m) 


COSY, TOCSY, J-res, 
HSQC 


16 


Rhamnose 


txCiH, aCjH, oiCjH, s(C4H, oiCsH, 
^CjH, /JCjH, /JC4H, /iCsH 


5.18 (d, 1.7 Hz), 3.91 (m), 3.79 (m), 3.43 (m), 3.85 (m), 3.92 (m), 
3.59 (m), 3.37 (m), 3.43 (m) 


COSY, TOCSY, J-res, 
HSQC 


17 


Rutinose 


Rha- C,H, Glc-C,H 


5.42 (d, 1.7 Hz), 4.50 (d, 7.4 Hz) 


COSY, TOCSY, J-res, 
HSQC 


18 


Rutinulose 


Rha- C,H, Fru-C,H 


5.25 (d, 1.7 Hz), 4.45 (d, 7.4 Hz) 


COSY, TOCSY, J-res, 
HSQC 


19 


Salvianic acid 


H-2, H-5, H-6, H2-7, H-8 


6.81 (d, 1.5 Hz), 6.86 (d, 8.0 Hz), 6.74 (dd, 8.0, 1.5 Hz), 2.85 (m), 
3.04 (m), 5.42 (m) 


COSY, TOCSY, J-res, 
HSQC 


20 


Salvianolic acid B 


H-2, H-3, H-5, H-6, H-/i, H-a, 

Ul i;" U A" u 7" Ul Q" Ul T" 
n-D , n-D , n2"/ , n-o , n-z , 

H-5'", H-6'", H-7'", H-8'" 


6.01 (d, 5.6 Hz), 4.33 (d, 5.6 Hz), 7.06, 6.92, 6.99 (d, 16.0 Hz), 
D.oD \0, 1 D.u nz), D.yD, o./D, z.yu imj, d.iu tmj, d.uu 
6.24 (d, 2.6), 6.42 (d, 8.0), 6.12 (dd, 8.0, 2.6 Hz), 2.98 (m), 
2.56 (m), 4.90 (m) 


COSY, TOCSY, J-res, 
nDtj^- 


21 


Rosmarinic acid 


H-2, H-5, H-6, H-ft H-o:, H-2', 
no , n o , ri2 ' 


6.98 (d, 1 .5 Hz), 6.72 (d,8.0 Hz), 6.89 (dd, 8.0, 1 .5 Hz), 

7 AO Iri 1ftnl-l-7l ft^lM lAnl-l?! ft Qn M 1 1-1-7 1 

/.'fy ^u, 1 D.u riz), D.DZ \o, lo.u nzj, o.yu \o, 1 . D nzj, 
6.64 (d, 8.0 Hz), 6.55 (dd, 8.0, 1.5 Hz), 2.91 (m), 3.13 (m) 


COSY, TOCSY, J-res, 
HSQC 


22 


Lithospermic acid 


uc Ufi U/-J UIT Liir' 

n-D, n-o, n-p, n-o!, n-z , n-D , 
H-6', H2-7', H-8', H-2", H-5", H-6" 


ft CTT lr\ Q K \-i-7\ ft 1^7 M Q c: I_I-t\ 7 T) (A 1 ft n I-I-7I ft lA 1 ft n 

D.DZ (u, o.D nz), o.D/ (u, o.D nZ), / .zz \u, iD.u nZj, o.dz \o, 10. u nz), 
6.74 (d, 1.5 Hz), 6.65 (d, 8.0 Hz), 6.52 (dd, 8.0, 1.5 Hz), 2.81 (rr), 
2.87 (m), 4.85 (m), 6.66 (d, 1.5 Hz), 6.55 (d, 8.0 Hz), 6.46 (dd, 8.0, 
1.5 Hz) 


{~r\QV Tnrcv / roc 
i-Udy, iuldy, J-res, 

HSQC 


23 


Salvianolic acid A 


H-5, H-6, H-p, H-a, H-2', H-5', 
H-6', H2-7', H-8', H-2", H-5", 
H-6", H-7", H-8" 


6.71 (d, 8.0 Hz), 7.05 (d, 8.0 Hz), 7.86 (d, 16.0 Hz), 6.20 (d, 16.0 Hz), 
6.66 (d, 1.5 Hz), 6.59 (d, 8.0 Hz), 6.48 (dd, 8.0, 1.5 Hz), 2.90 (m), 
3.01 (m), 5.1 1 (dd, 8.5, 3.5 Hz), 7.00 (d, 1.0 Hz), 6.69 (d, 8.0 Hz), 
o.oz (dd, o.o, l.D Hzj, o.bu (d, lo.o Hz), /.08 (d, ib.o Hz) 


COSY, TOCSY, J-res, 
HSQC 


24 


Procatechuic acid 


H-2, H-5, H-6 


7.41 {d, 2.0 Hz), 6.80 (d, 8.0 Hz), 7.45 (dd, 8.0, 2.0 Hz) 


COSY, TOCSY, J-res, 
HSQC 


25 


Procatechuic aldehyde 


H-2, H-5, H-6, CHO 


7.10 (d, 2.0 Hz), 6.80 (d, 8.0 Hz), 7.08 (dd, 8.0, 2.0 Hz), 9.63 (s) 


COSY, TOCSY, J-res, 
HSQC 


26 


4-Hydroxybenzoic acid 


H-3, 5, H-2, 6 


7.83 (d, 8.8 Hz), 6.91 (d, 8.8 Hz) 


COSY, TOCSY, J-res, 
HSQC 


27 


4-Hydroxycinnamic acid 


H-3, 5, H-2, 6, H-P, H-a 


7.54 (d, 9.0 Hz), 6.89 (d, 9.0 Hz), 7.60 (d, 16.0 Hz), 6.37 (d, 16.0 Hz} 


COSY, TOCSY, J-res, 
HSQC 


28 


Uridine 


H-5, H-6, Xyl-CiH 


5.88 (d, 8.1 Hz), 7.88 (d, 8.1 Hz), 5.91 (d, 4.2 Hz) 


COSY, TOCSY, J-res, 



HSQC 
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Table 1. Cont. 





Peak 


Metabolites 


Assignment 


J'H (multiplicity") 


Assigned with 


29 


Formate 


CHO 


8.46 (s) 


J-res, HSQC 


30 


5-(Hydroxymethyl)-2- 
furaldehyde 


H-3, H-4, CHO 


7.23 (d, 3.5 Hz), 6.51 (d, 3.5 Hz), 9.45 (s) 


COSY, TOCSY, J-res, 
HSQC 



^Multiplicity: singlet (s), doublet (d), triplet (t), doublet of doublets (dd), quintet (q), multiplet (m). 
doi:l 0.1 371 /journal.pone.Ol 0541 2.t001 



(Beijing, China). Salvianic acid and procatechuic acid were 
obtained from Zhongxin Innova Laboratories (Tianjin, RP 
China). Succinate and malonate were obtained from Dr. 
Ehrenstorfer GmbH (Augsburg, Germany). Fructose and glucose 
were purchased from Sigma (Aldrich, America). Rutinose was 
purchased from Hazard Communication (Tokyo, Japan). The 
purities of the compounds were all above 98%, using NMR 
analysis. Deuterium oxide (D2O, 99.9%) and sodium 3-trimethyl- 
silyl [2,2,3,3-6(4] propionate (TSP) were purchased from Cam- 
bridge Isotope Laboratories (Miami, FL, USA). D2O was used as 
internal lock; TSP was as internal standard for chemical shift 
calibration and quantification. 

Ethics 

No specific permission was required for the described field 
studies. The field locations are neither privately owned nor 
protected, and neither endangered nor protected species were 
involved. 

Sample preparation 

All the injection samples were subjected to freeze-drying. The 
dried powders (18 mg) were accurately weighed and dissolved with 
600 |iL of D,0 containing 0.58 mM TSP. Exactly 500 |iL of 
sample solution was transferred into a standard 5 mm NMR tube 
(Vineland, NJ, USA). No buffer was used due to the stable pH of 
injection. 

NMR measurements 

'H NMR spectra were acquired at 298 K on a Bruker AV III 
600 MHz NMR spectrometer (600.23 MHz for proton frequency) 
with a 5 mm broadband BBFO probehead. All pulse sequences 
were from Bruker pulse program library. A standard one 
dimensional composite pulse sequence (zgcppr) was employed to 
suppress the residual water signal. The 90° pulse width was 
adjusted to about 13 |J.m for each sample. Sixty-four scans were 
collected into 32k data points using a spectral width of 12335 Hz, 
a relaxation delay of 1.0 s and an acquisition time of 2.66 s. A 
0.3 Hz line-broadening function was applied to all spectra for 
Fourier transformation (FT) followed by phasing and baseline 
correction. 

For proton signals assignment purposes, a set of two dimen- 
sional (2D) spectra, including 'H-'H correlation spectroscopy 
(COSY), 'H-'H total correlation spectroscopy (TOCSY), 'H /- 
resolved (/-res), and 'H-'^'C heteronuclear single quantum 
coherence (HSQC), were acquired for selected samples and 
processed with similar parameters as described previously [35]. 

Spin-lattice relaxation time (Ti) values of the quantified protons 
of individual constituent and TSP were measured using a classical 
inversion recovery pulse sequence with 10 relaxation delays (t) 
ranging from 0.01 to 20 s. 



Quantification of the metabolites 

Because the intensity of a given 'H NMR signal is directiy 
proportional to its contributing number of protons, the amount of 
metabolites in DHI can be measured by the signal areas of given 
metabolites and an internal reference with known concentration. 
Thirteen metabolites in DHI were selected for quantification. The 
important parameters for data acquisition and processing of 'H 
NMR spectra must be set appropriately to obtain accurate and 
precise measurements [36]. Most of all, the relaxation delay T 
should be long enough to ensure complete relaxation for all the 
signals of interest. 'H NMR measurements are done in a longer 
acquisition time by choosing T^5x longest Ti. For shortening the 
acquisition time, H NMR spectra can be acquired in an 
incompletely relaxed condition, and the absolute concentrations 
should be calculated taking the T] values in consideration [37]. 
With the effective magnetization reading pulse of 90°, the 
quantification of chemical constituents in this study can be 
performed by using the following equation: 

_ Aj, ^-^^v{-7T[^')NrsP My ^^^^ 
^ Atsp \-exp(^-t/j^x'^ Mtsp ™^ 

where Px and Ptsp are the mass concentrations of metabolite and 
TSP, Ax and Aj-sp are the integral areas for targeted signal of 
metabolite and for methyl groups of TSP, Nx and Ntsp are the 
proton numbers of metabolite and of methyl groups of TSP, Mx 
and Mjsp are the molar masses of metabolite and TSP, and 
fTSP ^j^g spin-lattice relaxation times for proton X and methyl 
protons of TSP, respectively; t is total relaxation time (relaxation 
delay plus acquisition time). 

The quantitative 'H NMR method was checked for linearity, 
precision, repeatability, stability, and accuracy. Precision, repeat- 
ability and stability were calculated as the relative standard 
deviation (RSD). Recovery test was employed to determine the 
accuracy, and four typical compounds, alanine, glucose, salvianic 
acid and procatechuic aldehyde, were chosen to evaluate the 
average recovery. 

Assay for multivariate quality control 

36 batches of DHI samples were split into two phases. Phase I 
and Phase II. Phase I as training set included 25 batches of 
qualified products manufactured continuously in the years of 2012 
and 2013. Phase II as testing set included 6 batches of qualified 
products in 2012 and 2013, and 5 batches of expired products 
manufactured in 201 1. 

The quantitative resulting data of 1 3 metabolites was imported 
to R 3.0.2 software loaded with packages of MVA, MSQC and 
mixOmics (www.r-project.org) for multivariate statistical analysis. 
For simplifying the multivariate problem, principal component 
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Figure 2. Quantitative data of selected 13 metabolites. Boxplots of the quantitative data before and after centering and scaling in Phase I (A, 

B) and Phase II (C, D). 

doi:1 0.1 371 /journal.pone.01 0541 2.g002 



analysis (PCA) and independence principal component analysis 
(IPC A) were performed to reduce the dimensionality of data. The 
scores of principal components characterized the whole data were 
then imported to and hoteUing control charts to calculate 
the upper control limits (with 99% confidence ellipsoids). 
Moreover, as one of the requisites in control chart is the 
independence of the data, the independence of selected compo- 
nents in PCA and IPCA models were validated by autocorrelation 
fimction (ACF). The out-of-control samples in Phase II were 
examined by the upper control limits achieved from Phase I. 

Result and Discussion 

Proton signal assignments and chemical identification 

A representative 'H NMR spectrum of DHI was shown in 
Fig. 1 . The resonance signals were assigned to 30 metabolites based 
on the elucidation with extensive 2D NMR experiments ('H-'H 
COSY, 'H-'H TOCSY, 'H/-resolved, HSQC),the literature data 
in our former work [35], and in-house database. In the range of (5 
3.2-5.8, the spectrum is dominated by Smonosaccharides and 2 
disaccharides, including glucose, galactose, arabinose, fructose, 
rhamnose, rutinose, and rutinulose. In the high-field region [5 0.5- 
3.2), 8 amino acids (isoleucine, leucine, valine, threonine, alanine. 



proline, pyroglutamate, asparagine) and 3 organic acids (acetate, 
succinate, malonate) were observed. In the low-field region [S 5.8— 
10.0), 7 polyphenolic acids, including saivianic acid, salvianolic 
acid B, salvianolic acid A, rosmarinic acid, lithospermic acid, 
procatechuic acid and procatechuic aldehyde, together with- 
uridine and 5-(hydroxymethyl)-2-furaldehyde (5-HMF), were 
identified. Moreover, 3 organic acids (4-hydroxybenzoic acid, 4- 
hydroxycinnamic acid, and formate) were observed as well in the 
low-field region. The chemical shifts of the identified 30 
metabolites by'H NMR were listed in Table l.To our knowledge, 
7 saccharides and 6 organic acids were reported for the first time 
in DHI. Without the need of any sample pretreatment or pre- 
column derivatization, the established H NMR method provided 
an approach to determine 7 saccharides, 6 organic acids, 8 amino 
acids, 1 nucleoside, 1 carbohydrate derivatives (5-HMF) and 7 
polyphenolic acids in DHI simultaneously. 

Quantitative NMR analysis and method validation 

Due to the narrow chemical shift range of 'H NMR and 
frequent signal overlap, it is a challenge to quantify all the 
constituents in a mixture. In our study, 13 metabolites with fuUy 
separated signals were selected for quantification. In order to 
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Figure 3. Component contributions of PCA and IPCA. (A) Pareto chart of principal components (PCs) in PCA indicated the first three PCs in 
Phase I explained 80.99% variances. (B) Kurtosis measurements of all extracted independence principal components (IPCs) in IPCA with Phase I data. 
The kurtosis of IPCS was close to zero, so the first 7 components of IPCA were used to choose exact numbers of IPCs. (C) Kurtosis measurements of 
the first 7 IPCs showed that the kurtosis of IPCS was close to zero, and the first 2 components were sufficient with IPCA. 
doi:1 0.1 371 /journal.pone.01 0541 2.g003 



improve the efficiency of 'H NMR method, the spectra were 
acquired in an incompletely relaxed state. As a consequence, the 
spin-lattice relaxation time (Ti) value must be accurately measured 
and taken into account for quantitative analysis. The values 
were determined by the inversion-recovery experiments (Table 2). 
Accordingly, the absolute concentrations of the 13 metabolites 
were calculated from three parallel samples of each batch (Table 
SI in File SI). 

Linearity. 'H NMR as method itself is linear and no 
calibration is necessary for the determination of molar ratios of 
mixtures. Thus, the 13 metabolites in five dilferent molar ratios 
confirmed the linearity of NMR spectroscopy. Good linearity was 
achievable, as indicated by the equations and satisfactory 
correlation coefficients (r^). 

Precision. The intraday and interday precision was deter- 
mined by analyzing six replicates on the same day and on three 
consecutive days respectively. The intraday precision for the 
contents of 13 metabolites ranged from 0.20% to 0.89%, and the 
interday precision ranged from 0.29% to 1.49%. The RSD values 
were adequate and indicated the suitability of the method. 



Repeatability. Six samples prepared from the same batch 
showed RSD values ranging from 1.46% to 2.75%, indicating a 
high repeatability. 

Stability. One sample was analyzed to determine stability on 
three consecutive days. The RSD values of the analytes were in the 
range of 0.37% to 2.21%. 

Accuracy. Considering the limited volume of NMR tube and 
the cost of using deuterium reagents to dilute sample continuously, 
four metabolites of different types were employed, including 
alanine, glucose, salvianic acid and procatechuic aldehyde. The 
recovery was calculated as the ratio of the response of the selected 
four compounds in the spiked DHI samples against that of the 
standards at the same levels. The average recoveries were found to 
be 106.6% (±1.2), 106.6% (±2.5), 98.3% (±2.1) and 91.3% 
(±4.0) for alanine, glucose, salvianic acid and procatechuic 
aldehyde, respectively, indicating acceptable recovery. 

According to the results, the concentration of glucose was 
extremely high in DHI as shown in Fig. 2A and 2B, and the 
amount of saccharides, amino acids and organic acids represented 
about 60% of the total solid content of DHI (Table SI in File SI). 
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Figure 4. Multivariate quality control charts. The and Hotelling control charts of the first three principle components (PCs) in Phase I (A and 
B), and the first two independence principle components (IPCs) in Phase I (C and D) and Phase II (E and F). Red solid lines represented the upper 
control limits (UCL) of control charts with 99% confidence region. 
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Figure 5. Autocorrelation function (ACF) validation of PCA and IPCA. Correlograms for the first three principle components (PCI, PC2, PCS) 
and the first two independence principle components (IPC1, IPC2) in Phase I. Blue dashed lines represented the confidence bands of ACF. PCI fell 
outside of the confidence bands with autocorrelation or dependence, it was said that PCA model was not successful. 
doi:1 0.1 371 /journal.pone.01 0541 2.g005 



Control Charts based on PCA and IPCA 

With the enhancement iii quahty control of DHI, the analysis 
should be performed through a multivariate approach, that is, the 
above 13 metabolites must be analyzed together, not indepen- 
dently. The concentrations of 1 3 metabolites were mean centered 
and unit-variance scaled before being analyzed by PCA and IPCA 
(Fig. 2C and 2D). To avoid the loss of significant information, the 
percent specified of the principal components (PCs) cumulative 
proportion of explain variance is normally fixed on 80% in PCA 
model [29]. Thus, the first three PCs in Phase I (Fig. 3A) were 
selected to construct (Fig. 4A) and Hotelling T^ (Fig. 4B) 
control charts. Since the independence of the data is one of the 



requisites in control chart, we assessed the marginal independence 
of each necessary PC to indicate the model validation [38]. 
Correlograms (Fig. 5 PC1-PC3) showed that PCI fell outside of 
the confidence bands, which indicated that there was an evidence 
of autocorrelation or dependence of PC 1 , and PCA model using 
Phase I data was not achieved. 

In order to remove the autocorrelation effects of PCA, we 
employed IPCA to generate denoised and independent loading 
vectors [28]. The kurtosis measure of loading vectors was used to 
decide the number of independent principal components (IPCs). 
The kurtosis of all extracted IPCs was plotted in Fig. 3B, whereas 
the kurtosis of IPC7 was close to zero. By using the first 7 
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Figure 6. The confidence ellipses for quality control. Scatterplots for the first two IPC scores with the confidence ellipses in Phase 1 (A) and 
Phase II (B). The control ellipse used as process region was shown in green solid line, while the Hotelling control ellipse used as tolerance region 
was in red dashed line. The numbers represented the samples in Phase II listing in Table S2 in File SI, and the seventh sample fell outside the 
tolerance region. 
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Figure 7. The loading plots and values of IPCA. (A) The IPCA loading plots for the first two IPCs in Phase II, The numbers represented the 
samples in Phase II listing in Table S2 in File SI, The influence of metabolites was showing in red arrow lines, (B) Bar diagram showed the loading 
values of metabolites in first two independence principle components, 
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components of IPCA, the exactly choosing number of IPCs was 
obtained (Fig, 3C), Since tlie kurtosis of IPCS was close to zero, 
the first 2 components were sufficient with IPCA, The presence of 
autocorrelation was assessed as shown in Fig, 5 (IPC I and IPC2), 
which indicated that there was no evidence of relation between the 
adjacent observations. Therefore the original 13 dimension of our 
data had been reduced to a two-dimensional problem. 

Then the first two IPCs were taken to establish the in-control 
state (Phase I), According to (Fig, 4C) and Hotelling T^ 
(Fig, 4D) control charts with the 99% confidence region, the upper 
control limits (UCL) were determined as 9,21 and 11,56, 
respectively. The first two IPCs were consequently controlled 
through 2D ellipsoids (Fig, 6), The x"' control ellipse (UCL = 9,21) 
could be used as process region, and the T^ control ellipse with less 
restrictive (UCL= 11,56) was used as tolerance region. In both 
cases, all the points in Phase I fell inside the confidence ellipsoids 
(Fig, 6A), The samples in Phase II (Table S2 in File SI) were 
monitored by employing the UCLs of both and T^ charts 
obtained from Phase I (Fig, 4E and 4F), and the points of Phase II 
were added into 2D ellipsoids (Fig, 6B), The seventh sample fell 
outside the 99th confidence ellipsoids of both the process and 
tolerance regions, indicating the presence of out-of-control sample 
(batch 110402), The decomposition of T^ value showed that the 
out-of-control variability was associated to the IPC 1 since ^-value 
was equal to 0,0032 (Table 3). The same result was also obtained 
through IPCA loading plots as shown in Fig, 7, Although the 
loading plot could not exactly determine which metabolites were 
responsible for the variation, it still showed the contents of 

Table 3. Decomposition of the seventh sample in Phase II, 



succinate, malonate, glucose, fructose, salvianic acid and proto- 
catechuic aldehyde made more contributions to the independent 
loading vectors of IPC 1 , 

Conclusion 

Based on quantitative 'H NMR analysis, a reliable approach for 
simultaneous determination of amino acids, organic acids, 
saccharides, and botanic secondary metabolites of HMIs in one 
fingerprinting spectrum has been developed and validated by 
using Danhong injection as a model. The method had takenTi 
values into account when calculated the contents of feature 
metabolites, which allowed the assay with good linearity, precision, 
repeatability, stability and accuracy. Unlike HPLC fingerprinting 
methods, the H NMR approach has the significant advantages of 
less analysis time (about 5 min) without chromatographic separa- 
tion, and no requirement of standard materials used to quantita- 
tive analysis. In combination with IPCA, two kinds of multivariate 
control charts (x^ and HoteUing T^) were also successfully carried 
out for detecting off-test HMI samples by employing the 
independence principle components. The decomposition of T^ 
value and IPCA loading vectors can reflect the significance 
variations of overall metabolite profiling, although it cannot decide 
the mutative metabolites exactiy. The established multivariate 
models have the prospects for extracting sufficient characterization 
to monitor more feature metabolites in HMIs, 



Decomposition UCL p-Value 

[IJ 10.7286 7,5731 0,0032 

[2,1 0.9184 7.5731 0,3475 

[3,1 11.6469 11.5641 0,0000 

doi:l 0.1 371 /journal.pone.Ol 0541 2.t003 
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Supporting Information 

File SI Quantification of 13 chemical markers selected 
in the IH NMR spectra of DHI samples in Phase I (Table 
SI) and Phase II (Table S2). 

(DOC) 
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